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Recent developments in neutron star theory and observation are discussed. Based on 
modern nucleon-nucleon potentials more reliable equations of state for dense nuclear 
matter have been constructed. Furthermore, phase transitions such as pion, kaon and 
hyperon condensation, superfluidity and quark matter can occur in cores of neutron 
stars. Specifically, the nuclear to quark matter phase transition and its mixed phases 
with intriguing structures is treated. Rotating neutron stars with and without phase 
transitions are discussed and compared to observed masses, radii and glitches. The 
observations of possible heavy ~ 2Mq neutron stars in X-ray binaries and QPO's require 
relatively stiff equation of states and restrict strong phase transitions to occur at very 
high nuclear densities only. 



1 Introduction 

Neutron stars are complicated many-body systems of neutrons, protons, electrons 
and muons with possibly additional intricate phases of pion, kaon, hyperon or quark 
condensates. It is therefore appropriate that neutron stars are included at this "Xth 
Intl. Conf. on Recent Progress in Many-body Theories" . 

A brief history of the most important discoveries in this millennium concerning 
neutron stars listed in Table I. Most are well known except perhaps for the most 
recent ones discussed the following sections: Sec. 2 list neutron star masses from 
binary pulsars and X-ray binaries and specifically discuss recent masses and radii 
from quasi periodic oscillations in low mass X-ray pulsars. In Sec. 3 we turn 
to modern equation of states for neutron star matter with particular attention 
to the uncertainty in the stiffness of the equation of state at high densities and 
causality constraints. Sec. 4 attempts to give an up-to-date on possible phase 
transitions to kaon and pion condensates, hyperon and quark matter, superfluidity, 
etc. Sec. 5 contains the resulting structure of neutron stars, calculated masses 
and radii, and compares to observations. In Sec. 6 observational effects of glitches 
are described when phase transitions occur. Sec. 7 contain novel information on 
connections between supernova remnants and NO^ peaks in ice cores. X-ray bursts 
and thermonuclear explosions, gamma ray bursters and neutron star collapse to 
black holes. Finally, a summary and conclusion is given. For more details we refer 
toB. 



2 Observed neutron star masses and QPO's 

The best determined neutron star masses are found in binary pulsars and all lie in 
the range 1.35 ± O.O4M0 El. These masses have been accurately determined from 
variations in their radio pulses due to doppler shifts as well periastron advances of 
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Table 1: Chronological list of important developments related to neutron stars 



Year "Observers" Discovery 



1054 


Chinese 


record the Crab Supernova 


1572 


Tycho Brahe 


observes a Supernova 


1932 


Chadwick 


discovers the neutron 


1932 


Landau & others 


suggests the existence of neutron stars 


1934 


Baade & Zwicky 


connects supernovae to gravitational 
collapse ot stars to neutron stars 


1935 


Oppenheimer & VolkofF 


calculate first neutron star structures 


1946 


Gamow 


develops nucleosynthesis which 

requires heavier elements from Supernovae 


1967 


Bell & Hewish 


discover first pulsar 


1969 




pulsars in Crab and Vela Supernova remnants 


1973 


Hulse & Taylor 


discover first binary pulsar 


1987 


Neutrino detectors 


collect 19 neutrinos from SN-1987A 


1995 


Nijmegen data base 


compilation of ^5000 NN cross sections 
leads to "modern" NN potentials and EoS 


1996 


RXTE 


kHz oscillations (QPO) in X-ray binaries 


1997 


BeppoSAX 


Gamma Ray Burst with afterglow at z^l 
Beaming and syncroton radiation observed 


1998 


Supernova remnants 


match dating of polar ice core iVO;^-peaks 



their close elliptic orbits that are strongly affected by general relativistic effects. One 
ception is the nonrelativistic pulsar PSR J1012+5307 of mass^Af = (2.1±0.8)Mq 



Several Xpray binary masses have been measured of which the heaviest are 
Vela X-l withy M = (1.9 ± 0.2)Mq and Cygnus X-2 withQ M = (1.8 ± OA)Mq. 
Their Kepler orbits are determined by measuring doppler shift of both the X-ray 
binary and its companion. To complete the mass determination one needs the 
orbital inclination which is determined by eclipse durations, optical light curves, or 
polarization variationsEI. 

The recent discovery of high-frequency brightness oscillations in low-mass X- 
ray binaries movides a promising new method for determining masses and radii of 
neutron starsQ. The kilohertz quasi-periodic oscillations (QPO) occur in pairs and 
are most likely the orbital frequencies 



i^QPO^ il/2TT) ^GM/Rl,, (1) 

of accreting matter in Keplerian orbits around neutron stars of mass M and its 
beat frequency with the neutron star spin, i^Qpa^ v^. The accretion can for a few 
QPO's be tracked to its innermost stable orbit 0i3 

Ems = 6GM/c^ . (2) 



All uncertainties given here are 95% conf. limits or ~ 2a 
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For slowly rotating stars the resulting mass is from Eqs. 

kHz 

M~2.2Mq . (3) 

vqpo 

For example, the maximum frequency of 1060 Hz upper QPO observed in 4U 1820- 
30 gives M ~ 2.25M0 after correcting for the Vs — 275 Hz neutron star rotation 
frequency. If the maximum QPO frequencies of 4U 1608-52 {vqpo = 1125 Hz) 
and 4U 1636-536 {vqpo — 1228 Hz) also correspond to innermost stable orbits the 
corresponding masses are 2.IM0 and 1.9Mq. Evidence for the innermost stable 
orbit has been found for 4U 1820-30, where vqpo display a distinct saturation with 
accretion rate indicating that orbital frequency cannot exceed that of the innermost 
stable orbit. However, one caveat is that large accretion leads to radiation that can 
slow down the rotation of the accreting matter. More observations are needed before 
a firm conclusion can be made. 

Large neutron star masses of order ~ 2Mq would restrict the equation of state 
(EoS) severely for dense matter as addressed in the following. 



3 Modern nuclear equation of states 

Recent models for the nucleon-nucleon (NN) interaction, based on the compilation 
of more than 5000 NN cross sections in the Nijmegen data bank, have reduced the 
uncertainty in NN potentials. The last Indiana run at higher momenta will further 
reduce uncertainties in NN interactions. Including many-body effects, three-body 
forces, relativistic effects, etc., the nuclear EoS have been constructed with reduced 
uncertainty allowing for more reliable calculations of neutron star properties 
Likewise, recent realistic effective interactions for nuclear matter obeying causality 
at high densities, constrain-the EoS severely and thus also the maximum masses of 
neutron stars. We have in EJ elaborated on these analyses by incorporating causality 
smoothly in the EoS for nuclear matter allowing for first and second order phase 
transitions to, e.g., quark matter. 

For the discussion of the gross properties jOf neutron stars we will use the optimal 
EoS of Akmal, Pandharipande, & Ravenhallll3 (specifically the Argonne V18 + Sv+ 
UIX* model- hereafter APR98), which is bascd_on the most recent models for the 
nucleon-nucleon interaction, see Engvik et al. t3 for a discussion of these models, 
and with the inclusion of a parametrized three-body force and relativistic boost 
corrections. The EoS for nuclear matter is thus known to some accuracy for densities 
up to a few times nuclear saturation density no = 0.16 fm~^. We parametrize 
the APR98 EoS by a simple form for the compressional and symmetry energies 
that gives a good fit around nuclear saturation densities and smoothly incorporates 
causality at high densities such that the sound speed approaches the speed of light. 
This requires that the compressional part of the energy per nucleon is quadratic in 
nuclear density with a minimum at saturation but linear at high densities 

£ = E,o,np{n) + S{n){l - 2xf 
7; — 2 — s 

^£ou- + W(l-2a;)2. (4) 

1 + su 



3 




Here, n ^ Up + n„ is the total baryon density, x — Up/n the proton fraction and 
u = n/no is the ratio of the baryon density to nuclear saturation density. The 
compressional term is in Eq. (^) parametrized by a simple form which reproduces 
the saturation density and the binding energy per nucleon Eq — 15.8MeV at uq of 
APR98. The "softness!' parameter s ~ 0.2, which gave the best fit to the data of 
APR98 is determined El by fitting the energy per nucleon of APR98 up to densities 
of n '--^ 4no. For the symmetry energy term we obtain 5*0 = 32 MeV and 7 = 0.6 
for the best fit. The proton fraction is given by /^-equilibrium at a given density. 

The one unknown parameter s expresses the uncertainty in the EoS at high 
density and we shall vary this parameter within the allowed limits in the following 
with and without phase transitions to calculate mass, radius and density relations 
for neutron stars. The "softness" parameter s is related to the incompressibility of 
nuclear matter as K^= 18fo/(l + s) — 200MeV. It agrees with the poorly known 
experimental value E3, Kq ~ 180 — 250MeV which does not restrict it very well. 
From {vs/cY — dP/d{n£), where P is the pressure, and the EoS of Eq. (||), the 
causality condition Cs < c requires 

0.13, (5) 

where nin is the mass of the nucleon. With this condition we have a causal EoS 
that reproduces the data of APR98 at densities up to 0.6 ~ 0.7 fm~'^. In contrast, 
the EoS of APR98 becomes superluminal at n « 1.1 fm~^. For larger s values the 
EoS is softer which eventually leads to smaller maximum masses of neutron stars. 
The observed M ~ I.4M0 in binary pulsars restricts s to be less than 0.4 — 0.5 
depending on rotation as shown in calculations of neutron stars below. 

In Fig. |l| we plot the sound speed (ws/c)^ for various values of s and that result- 
ing from the microscopic calculation of APR98 for /3-stable pn-matter. The form 
of Eq. (^, with the inclusion of the parameter s, provides a smooth extrapolation 
from small to large densities such that the sound speed Vg approaches the speed of 
light. For s = 0.0 (s = 0.1) the EoS becomes superluminal at densities of the order 
of 1 (6) fm-3. 

The sound speed of Kalogera & Baym is also plotted in Fig. |^. It jumps 
discontinuously to the speed of light at a chosen density. With this prescription 
they were able to obtain an optimum upper bound for neutron star masses and 
obey causality. This prescription was also employed by APR98. The EoS is thus 
discontinuously stiffened by taking = c at densities above a certain value ric 
which, however, is lower than = 5no where their nuclear EoS becomes super- 
luminal. This approach stiffens the nuclear EoS for densities ric < n < Ug but 
softens it at higher densities. Their resulting maximum masses lie in the range 
2.2Mq^M^2.9Mq. Our approach however, incorporates causality by reducing 
the sound speed smoothly towards the speed of light at high densities. Therefore 
our approach will not yield an absolute upper bound on the maximum mass of a 
neutron star but gives reasonable estimates based on modern EoS around nuclear 
matter densities, causality constraints at high densities and a smooth extrapolation 
between these two limits (see Fig. |l|). 

At very high densities particles are expected to be relativistic and the sound 
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Figure 1: (vs/cY for /3-stable pn-matter for s = 0, 0.1, 0.2, 0.3, the results of APR98, and for the 
patched EoS of Kalogera &; Baym which shows a discontinuous (ua/c)^. 



speed should be smaller than the speed of light, v1 ~ Consequently, the EoS 

should be even softer at high densities and the maximum masses we obtain with 
the EoS of (^) are likely to be too high estimates. 

4 Phase transitions 

The physical state of matter in the interiors of neutron stars at densities above a 
few times normal nuclear matter densities is essentially unknown and many first 
and second order phase transitions have been speculated upon. 

Kaon condensation 

Kaon condensation in dense matter was suggestfid by Kaplan and Nelson 111, and 
has been discussed in many recent publications EJIla. Due to the attraction between 
K~ and nucleons its energy decreases with increasing density, and eventually if it 
drops below the electron chemical potential in neutron star matter in /3-equilibrium, 
a Bose condensate of will appear. It is found that A'^'s condense at densities 
above ~ 3 — 4po, where po = 0.16 fm~^ is normal nuclear matter density. This is 
to be compared to the central density of ^ ApQ for a neutmn star of mass IAMq 
according to the estimates of Wiringa, Fiks and Fabrocini Ej using realistic models 
of nuclear forces. 

In neutron matter at low densities, when the interparticle spacing is much larger 
than the range of the interaction, rp S> i?, the kaon interacts strongly many times 
with the same nucleon before it encounters and interacts with another nucleon. 
Thus one can use the scattering length as the "effective" kaon-neutron interaction, 
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a^-N — — 0.41fni, where we ignore the minor proton fraction in nuclear matter. 
The kaon energy deviates from its rest mass by the Lenz potential 

27r . , 

rriR 

which is the optical potential obtained in the impulse approximation. If hadron 
massesJurthermore decrease with density the condensation will occur at lower den- 
sities M. 

At high densities when the interparticle spacing is much less than the range of 
the interaction, tq <C R, the kaon will interact with many nucleons on a distance 
scale much less than the range of the interaction. The kaon thus experiences the 
field from many nucleons and the kaon energy deviates from its rest mass by the 
Hartree potential: 

UJHartree = TTlK + TlNM J VK-N{r)d^r , (7) 

As shown in Ref. 0, the Hartree potential is considerably less attractive than the 
Lenz potential. Already at rather low densities, when the interparticle distance 
is comparable to the range of the KN interaction, the kaon-nucleon and .niicleon- 
nucleon correlations conspire to reduce the K~ N attraction significantly 113113. This 
is also evident from Fig. |^ where the transition from the low density Lenz poten- 
tial to the high density Hartree potential is calculated by solving the Klein-Gordon 
equation for kaons in neutron matter in the Wigner-Seitz cell approximation. Re- 
sults are for square well if ^ A^-potentials of various ranges R. For the measured 
K~n scattering lengths and reasonable ranges of interactions the attraction is re- 
duced by about a factor of 2-3 in cores of neutron stars. Relativistic effects further 
reduce the attraction at high densities. Consequently, a kaon condensate is less 
likely in neutron stars due to nuclear correlations. 

If the kaon condensate occurs a mixed phase of kaon condensates and ordinary 
nuclear matter may coexist in a mixed phase l3 depending on surface and Coulomb 
energies involved. The structures would be much like the quark and nuclear matter 
mixed phases described above. 



J;.. 2 Pion condensation 

Pion condcnaation is like kaon condensation possible in dense neutron star matter 
(see, e.g., l^'lII). If we first neglect the effect of strong correlations of pions with the 
matter in modifying the pion self-energy, one finds it is favorable for a neutron on 
the top of the Fermi sea to turn into a proton and a tt^ when 

^J■n - fJ-p ^ ^J.e > m^- , (8) 

where m^- — 139.6 MeV is the rest mass. As discussed in the previous subsec- 
tion, at nuclear matter saturation density the electron chemical potential is ^ 100 
MeV and one might therefore expect the appearance of 7r~ at a slightly higher den- 
sity. One can however not neglect the interaction of the pion with the background 
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Figure 2: Kaon energy as function of neutron density is shown by full curves for various ranges of 
the K~n potentials -R = 0.4,0.7, 1.0 fm. Nuclear correlations change the low density Lenz result 
(Eq. (^, dotted curve) to the high density Hartree result (Eq.(|^), dashed curves). The electron 
chemical potential /ie from the EoS of Eq. (U) with 5 = 0.2 is shown with and without (lower and 
upper dotted curves) a transition to a mixed phase of quark matter for B = 100 MeVfm"''. 

matter. Such interactions can enhance the pion self-energy and thereby tha,pion 
threshold density, and depending on the chosen parameters, see again Ref. c3, the 
critical density for pion condensation may vary from hq to 4no. These matters are 
however not yet settled in a satisfying way, and models with strong nucleon-nucleon 
correlations tend to suppress both the ttNN and ttAN interaction vertices so that 
a pion condensation in neutron star matter does not occur. 

A 7r° condensate may also form as recently suggested by Akmal et al. E3 and 
appears at a density of ^ 0.2 fm~'^ for pure neutron matter when the three-body 
interaction is included, whereas without Vijk it appears at much higher densities, 
i.e. ^ 0.5 fm~^. The tt'~' are virtual in the same way as the photons keeping a solid 
together are virtual. 



4-3 Hyperon matter 

Condensates of hyperons A, 1]^'"'+, ... and Deltas A~''''+'++ also appear when their 

chemical potential exceeds their effective mass in matter. In /3-equilibrium the 
chemical potentials are related by 

A^ri + A*e = A^s- = A^A- (9) 

M« = /"A = MSO = MAO (10) 

^^p = A*A+ • (11) 
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The appears via weak strangeness non-conserving interactions +ri — > +Ve, 
when /Lis- > co-^- and A hyperons when /iA > If we neglect interactions, one 
would expect the S~ to appear at lower densities than the A, even though E~ is 
more massive due to the large electron chemical potential. The threshold densities 
for noninteracting E~,A~,A are relatively low indicating that these condensates 
are present in cores of neutron stars if their interactions can be ignored. 

Hyperon energies are, however, strongly affected in dense nuclear matter by 
interactions and correlations with the nucleons and other hyperons in a condensate. 
If hyperons have the short range repulsion and three-body interactions as nucleons, 
a condensate of hyperons become less likely. 

4-4 Quark matter 

Eventually at high densities, we expect hadrons to deconfine to quark matter or, in 
other words, chiral symmetry to be restored. This transition has been extensively 
studied by the Bag model equation of state (EoS) for quark matter which leads 
to a first order phase transition from hadronic to quark matter at a density and 
temperature determined by the parameters going into both EoS. In the bag model 
the quarks are assumed to be confined to a finite region of space, the so-called 'bag', 
by a vacuum pressure B. Adding the Fermi pressure and interactions computed to 
order as = /An, where g is the QCD coupling constant, the total pressure for 
three massless quarks of flavor f ^ u,d, s, is 

P^-P^{l^-as)-B + P„ + P^, (12) 

where Pg^^ are the electron and muon pressure, e.g., = /ig/127r^. A Fermi gas of 
quarks of flavor i has density Ui = kp^/Tr"^, due to the three color states. A finite 
strange quark mass have minor effect on the EoS since quark chemical potentials 
Hq^ruM/'i typically are much larger. The value of the bag constant B is poorly 
known, and we present results using two representative values, B = 150 MeVfm"'^ 
and B — 200 MeVfm"^. We take as — 0.4. However, similar results can be obtained 
with smaller as and larger B. 

The quark and nuclear matter mixed phase has continuous pressures and den- 
sities due to the general Gibbs criteria for two-component systems EZl. There are 
no first order phase transitions but at most two second order phase transitions. 
Namely, at a lower density, where quark matter first appears in nuclear matter, 
and at a very high density (if gravitationally stable), where all nucleons arc finally 
dissolved into quark matter. This mixed phase does, however, not include local 
surface and Coulomb energies of the quark and nuclear matter structures. If the 
interface tension between quark and nuclear matter is too large, the mixed phase 
is not favored energetically" due to surface and Coulomb energies associated with 
forming these structures o. The neutron star will then have a core of pure quark 
matter with a mantle of nuclear matter surrounding it and the two phases are co- 
existing by a first order phase transition or Maxwell construction. For a small or 
moderate interface tension the quarks are confined in droplet, rod- and plate-like 
structures as found in the inner crust of neutron stars (see t3 and Fig. ||). 
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Figure 3: The quark and nuclear matter structure in a quarter of a typical IAMq solar mass 
neutron star. The typical sizes of structures are a few Fermi's but have been scaled up by about 

16 orders of magnitudes to be seen. 

4-5 Superfluidity in baryonic matter 

The presence of neutron superfluidity in the crust and the inner part of neutron stars 
are considered weh estabhshed in the physics of these compact stellar objects. In 
the low density outer part of a neutron star, the neutron superfluidity is expected 
mainly in the attractive ^5*0 channel. At higher density, the nuclei in the crust 
dissolve, and one expects a region consisting of a quantum liquid of neutrons and 
protons in beta equilibrium. The proton contaminant should be superfluid in the 
^5*0 channel, while neutron superfluidity is expected to occur mainly in the coupled 
^P2-^F2 two-neutron channel. In the core of the star any superfluid phase should 
finally disappear. 

Dilute Fermi systems can now be studied as atomic gases recently have been 
cooled dopim to nanokelvin similar to Bose-Einstein condensates. Degeneracy was 
observed Ej and BGS, gaps are currently searched for. According to Gorkov and 
Melik-Bharkuderov t3 the gap is 



in the dilute limit where the Fermi momentum times the scattering length is small, 
kF\a\ < 1. 

Recently color superconductivity in quark matter has been taken up again since 
the quark color interaction has been conjectured to be large. Correspondingly large 
gaps of order tens of MeV are found. If the strange quark mass is sufficiently small 
color-flavor locking occur between up and down or strange quarks. We refer to T. 
Schaefer these proceedings for further details. 




7/3 



(13) 
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5 Calculated neutron star masses and radii 



In order to obtain the mass and radius of a neutron star, we have solved the Tolman- 
Oppenheimer-Volkov equation with and without rotational corrections. The equa- 
tions of state employed are given by the pn-matter EoS with s = 0.13, 0.2, 0.3, 0.4 
with nucleonic degrees of freedom only. In addition we have selected two represen- 
tative values for the bag-model parameter B, namely 150 and 200 MeVfm^'^ for our 
discussion on eventual phase transitions. The quark phase is linked with our pn- 
matter EoS from Eq. (|j)-.with s = 0.2 through cither a mixed phase construction or 
a Maxwell constructional. For B — 150 MeVfm^'^, the mixed phase begins at 0.51 
fm~^ and the pure quark matter phase begins at 1.89 fm~^. Finally, for B = 200 
MeVfm"^, the mixed phase starts at 0.72 fm"'^ while the pure quark phase starts at 
2.11 fm^'^. In case of a Maxwell construction, in order to link the pn and the quark 
matter EoS, we obtain for B = 150 MeVfm"'^ that the pure pn phase ends at 0.92 
fm~'^ and that the pure quark phase starts at 1.215 fm^'^, while the corresponding 
numbers for B = 200 MeVfm^^ are 1.04 and 1.57 fm^^. 

None of the equations of state from either the pure pn phase or with a mixed 
phase or Maxwell construction with quark degrees of freedom, result in stable con- 
figurations for densities above ^ IOtiq, implying thereby that none of the stars have 
cores with a pure quark phase. The EoS with pn degrees of freedom have masses 
M^2.2M0 when rotational corrections are accounted for. With the inclusion of the 
mixed phase, the total mass is reduced since the EoS is softer. However, there is 
the possibility of making very heavy quark stars for very small bag constants. For 
pure quark stars_there is only one energy scale namely B which provides a homology 
transformation cB and the maximum mass is Mmax = 2.OM0(58MeVfm^^/i?)^/^ 
(for as = 0). However, for _B^58MeVfm^'^ a nuclear matter mantle has to be added 
and for B <,58MeVfm~'^ quark matter has lower energy per baryon than ^^Fe and is 
thus the ground state of strongly interacting matter. Unless the latter is the case, 
we can thus exclude the existence of such 2.2 — 2.3M0 quark stars. 

In Fig. H we show the mass-radius relations for the various equations of state. 
The shaded area represents the allowed masses and radii for i'qpo = 1060 Hz of 4U 
1820-30. Generally, 

where the lower limit insures that the star is not a black hole, and the upper limit 
that the accreting matter orbits outside the star, R < Rort- Furthermore, for the 
matter to be outside the innermost stable orbit, R > Rms = GGM, requires that 

M < ^ + 0^ 2.2Mo(l + 0.75j)^ , (15) 

12\/6ttGi^qpo vqpo 

where j = 2'kcVsI /JVP is a dimensionless measure of the angular momentum of 
the star with moment of inertia /. The upper limit in Eq. (^) is the mass when 
vqpo corresponds to Ahe innermost stable orbit. This is the case for 4U 1820- 
30 since according to vqpo saturates at ~ 1060 Hz with increasing count rate. 
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Figure 4: Neutron star masses vs. radius for the EoS of Eq. (|4|) with softness 8=0.13,0.2,0.3,0.4, 
with increasing values of s from top to bottom for the full curves. Phase transitions decrease the 
maximum mass whereas rotation increases it. The shaded area represents the neutron star radii 
and masses allowed (see text and Eqs. 1-3) for orbital QPO frequencies 1060 Hz of 4U 1820-30. 

The corresponding neutron star mass is M ~ 2.2 — 2.3AfQ which leads to several 
interesting conclusions as seen in Fig. ^. Firstly, the stifFest EoS allowed by causality 
(s ~ 0.13 — 0.2) is needed. Secondly, rotation must be included which increase the 
maximum mass and corresponding radii by 10-15% for ~ 300 Hz. Thirdly, a 
phase transition to quark matter below densities of order ~ 5no can be excluded, 
corresponding to a restriction on the bag constant 5^200 MeVfm^"^. 

These maximum masses are smaller than those of APR98 and Kalogera & 
Baym who, as discussed above, obtain upper bounds on the mass of neutron stars 
by discontinuously setting the sound speed to equal the speed of light above a 
certain density, ric- By varying the density ric = 2 — > Stiq the maximum mass drops 
from 2.9 — *■ 2.2M0. In our case, incorporating causality smoothly by introducing 
the parameter s in Eq. (^, the EoS is softened at higher densities in order to obey 
causality, and yields a maximum mass lower than the 2.2M0 derived in APR98 for 
nonrotating stars. 

If the QPOs are not from the innermost stable orbits and one finds that even 
accreting neutron stars have small masses, say like the binary pulsars, M^IAMq, 
this may indicate that heavier neutron stars are not stable. Therefore, the EoS 
is soft at high densities or a phase transition occurs at a few times nuclear 

matter densities. 
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6 Glitches and phase transitions in rotating neutron stars 

Younger pulsars rotate and slow down rapidly. Some display sudden speed ups 
referred to as glitches. The glitches observed in the Crab, Vela, and a few other 
pulsars are probably due to quakes occurring in solid structures^uch as the crust, 
superfluid vortices or possibly the quark matter lattice in the core £3. As the rotating 
neutron star gradually slows down and becomes less deformed, the rigid component 
is strained and eventually cracks/quakes and changes its structure towards being 
more spherical. 

The moment of inertia of the rigid component, /c, decreases abruptly and its 
rotation and pulsar frequency increases due to angular momentum conservation 
resulting in a glitch. The observed glitches are very small Af2/il ~ 10~^. The 
two components slowly relaxate to a common rotational frequency on a time scale 
of days (healing time) due to superfluidity of the other component (the neutron 
liquid). The healing parameter Q = Ic/hot measured in glitches reveals that for the 
Vela and Crab pulsar about ~3% and ^96% of the moment of inertia is in the rigid 
component respectively. 

If the crust were the only rigid component the Vela neutron star should be 
almost all crust. This would require that the Vela is a very light neutron star - 
much smaller than the observed ones which all are compatible with ~ IAMq. If 
we by the lattice component include not only the solid crust but also the protons in 
nuclear matter (NM) (which is locked to the crust due to magnetic fields), superfluid 
vortices pinned to the crust and the solid QM mixed phase 

Ic — Icrust ^ Ip+ Isv + IqM , (16) 

we can better explain the large Ic for the Crab. The moment of inertia of the 
|taixed phase is sensitive to the EoS's used. For example, for a quadratic NM EoS 
Ea decreasing the Bag constant from 110 to 95 MeVfm""^ increases Ic/Itotai from 
~ 20% to ^ 70% for a 1.4Mq neutron star - not including possible vortex pinning. 
The structures in the mixed phase would exhibit anisotropic elastic properties, being 
rigid to some shear strains but not others in much the same way as liquid crystals. 
Therefore the whole mixed phase might not be rigid. 

As the neutron star slows down, pressures and densities increase in the center 
and a first order phase transition may occur. Does it leave any detectable sisnal 
at the corresponding critical angular velocity f2o ? As described in detail in E3'EJ 
the general relativistic equations for slowly rotating stars can be solved even with 
first order phase transitions since only the monopole is important. The resulting 
moment of inertia have the characteristic behavior for fi^irio 

/ = /o (l + cifi^ - C2{nl - n'f'^ + . (17) 

Here, Ci and C2 are small parameters proportional to the density difference between 
the two phases; however, C2 = for > JIq. 

In order to make contact with observation, the temporal behavior of angular 
velocities must be considered. The pulsars slow down at a rate given by the loss 
of rotational energy which one usually assumes is proportional to the rotational 
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angular velocity to some power (for dipole radiation rt = 3) 



I (im^) ^ (18) 

With the moment of inertia given by Eq. ( [l^ ) the decreasing angular velocity 
can be found. The corresponding braking index depends on the second derivative 
I" — dI/(Pil of the moment of inertia and thus diverges as H. approaches SIq from 
below 

n{n) = ^^n~cin^ + C2^^—. (19) 

The observational braking index n{Q) exhibits a characteristic behavior different 
from the theoretical braking index n in case of a first order phase transition. 

The critical angular velocities depend on the EoS and the critical densities. For 
best detection one would want the transition to occur for rapidly rotating pulsars 
such as millisecond pulsars, X-ray binaries or young neutron stars only a few years 
or centuries old. As pulsars slow down over a million years, their central densities 
span a wide range of order several tiq. As we are interested in time scales of years, 
we must instead study the ^ 1000 pulsars available. By studying the corresponding 
range of angular velocities for the sample of different star masses, the chance for 
encountering a critical angular velocity increases. Eventually, one may be able to 
cover the full range of central densities and find all first order phase transitions up 
to a certain size determined by the experimental resolution. 



7 Other recent surprises related to neutron stars 

7.1 Connecting NO^ peaks in ice core samples to supernovae 

A curious connection between four NO^ peaks from south pole ice core samples 
and supernovae has recently been made. A_small nearby but well hidden supernova 
remnant RX J0852. 0-4622 was discovered l3 in the southeast corner of the older Vela 
supernova remnant. Estimates of its distance and age seem toiie compatible with 
the 4th NO^ peak in the 20 year old south pole core samples l3. The other three 
already agreed with the historical Crab, Tycho and Kepler supernovae within the ice 
dating accuracy of ±20 years, a millennium back. This seems to indicate that nearby 
supernova explosions can affect the climate on earth and lead to geophysical signals. 
The radiation does produce NO^ in the upper atmosphere but the quantitative 
amount of fall out on the poles by northern/southern light cannot be estimated. 
NO^ peaks were not found in the Vostok og Greenland drills possibly because the 
aurorae do not have fallout on these sites. Further drillings are needed to confirm 
the connection - in particular deeper ones which should also include the 1006 and 
1054 supernovae. Mass extinctions may be caused by nearby supernova explosions 
after all. 
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1.2 X-ray bursts and thermonuclear explosions 

Slow accretion from a small mass companion ^2Mq generates a continuous back- 
ground of X-rays. Each nucleon radiates its gravitational energy of ~ m„GA//i? ~ 
100 MeV. After accumulating hydrogen on the surface, pressures and temperatures 
become sufficient to trigger an irregular runaway thermonuclear explosion every few 
hours or so seen as an X-ray burst (see, e.g., a for a review). The energy involved 
is typical nuclear binding energies ~ 1 MeV, i.e., a percent of the time integrated 
background. 

The microsecond time resolution in these X-ray spectra allow for Fourier trans- 
forms in millisecond time intervals which is much shorter than the burst duration 
of order a few seconds. In the case of 4U 1728-34 the power analysis shows a peak 
at the neutron star spin frequency at 364,Hz which, however, decreases to 362 Hz 
during the first 1-2 seconds of the burst □. A simple explanation is that the ther- 
monuclear explosion elevates the surface of the neutron star. Conserving angular 
momentum, L cx MR^v, leads to a decrease in rotation by 

— --2—. 20 
v R 

With a frequency change of Aiy ~ — 2 Hz and typical neutron star radii of order R ~ 
10 km, we find an elevation of order AR ~ 20 m, which is roughly in agreement with 
expectations but much less than on earth due to the much stronger gravitational 
fields on neutron stars. 



7.3 Gamma Ray Bursters and hypernovae 

The recent discovery of afterglow in Gamma Ray Bursters (GRB) allows determi- 
nation of the very high redshifts {z > I) and thus the enormous distance and energy 
output E ~ 10^^ ergs in GRB iOsotropically emitted. Very recently evidence for 
beaming or jets has been found t3 corresponding to "only" E ~ lO^'^ ergs. Candi- 
dates for such violent events include neutron star mergers or a special class of type 
Ic supernova (hypernovae) where cores collapse to black holes. Indication of such 
connections are brought by recent observations of a bright supernova coinciding 
with GRB 980326. Binary pulsars are rapidly spiraling inwards and will eventually 
merge and create a gigantic explosion and perhaps collapse to a black hole. From 
the number of binary pulsars and spiral rates, one estimates about one merger per 
million year per galaxy. With ~ 10^ galaxies at cosmological distances z^l this 
gives about the observed rate of GRB. However, detailed calculations have problems 
with baryon contamination, i.e. the baryons ejected absorb photons in the relativis- 
tically expanding photosphere. Accreting black holes have also been suggested to 
act as beaming GRB. 

So far, the physics producing these GRB is not understood. The time scales 
and the enormous power output points towards neutron star or black hole objects. 
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8 Summary 



Modern nucleon-nucleon potentials have reduced the uncertainties in the calculated 
EoS. Using the most recent realistic effective interactions for nuclear matter of 
APR98 with a smooth extrapolation to high densities including causality, the EoS 
could be constrained by a "softness" parameter s which parametrizes the unknown 
stiffness of the EoS at high densities. Maximum masses have subsequently been 
calculated for rotating neutron stars with and without first and second order phase 
transitions to, e.g., quark matter at high densities. The calculated bounds for 
maximum masses leaves two natural options when compared to the observed neutron 
star masses: 

• Case I: The large masses of the neutron stars in QPO 4U 1820-30 (M = 
2.?,Mq), PSR J1012+5307 (M = 2.1±GAMq), Vela X-1 (M = I.Q±{).IMq), 
and Cygnus X-2 (M = 1.8 ± Q.2Mq), are confirmed and complemented by 
other neutron stars with masses around ~ 2Mq. The EoS of dense nuclear 
matter is then severely constrained and only the stiffest EoS consistent with 
causality are allowed, i.e., softness parameter 0.13 < s^O.2. Furthermore, any 
significant phase transition at densities below < 5no can be excluded. 

That the radio binary pulsars all have masses around 1.4Mq is then probably 
due to the formation mechanism in supernovae where the Chandrasekhar mass 
for iron cores are ~ 1.5Mq. Neutron stars in binaries can subsequently acquire 
larger masses by accretion as X-ray binaries. 

• Case II: The heavy neutron stars prove erroneous by more detailed observa- 
tions and only masses like those of binary pulsars are found. If accretion does 
not produce neutron stars heavier than ^IAMq, this indicates that heavier 
neutron stars simply are not stable which in turn implies a soft EoS, either 
s > 0.4 or a significant phase transition must occur already at a few times 
nuclear saturation densities. 

Surface temperatures can be estimated from spectra and from the measured 
fluxes and known distances, one can extract the surface area of the emitting spot. 
This gives unfortunately only a lower limit on the neutron star size, R. If it becomes 
possible to measure both mass and radii of neutron stars, one can plot an obser- 
vational (M, R) curve in Fig. (Q) , which uniquely determines the EoS for strongly 
interacting matter at low temperature. 

Pulsar rotation frequencies and glitches are other promising signals that could 
reveal phase transitions. Besides the standard glitches also giant glitches were 
mentioned and in particular the characteristic behavior of angular velocities when 
a first order phase transition occurs right in the center of the star. 

It is impossible to cover all the interesting and recent developments concerning 
neutron stars in these proceedings so for more details we refer to El and Refs. therein. 
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